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It is shown that the rate limiting time constant for the for- 
mation of the Kondo state in a quantum dot can be extracted 
analytically from known perturbation theoretic results. The 
prediction obtained is verified via numerical simulations in 
the noncrossing approximation. 



The basic time scale for equilibration processes in a 
Kondo system has eluded researchers for more than thirty 
five years. The Kondo effect has recently been unambigu- 
ously observed in quantum dots [0J^]. These are proto- 
type Kondo systems, which are also important in their 
own right, and appropriate time dependent experiments 
should soon appear. The observation of Kondo resonance 
in the tunneling from an STM tip via a surface adsorbate 
atom has further increased the excitement. Impu- 
rity models incorporating Kondo physics form one of the 
basic building blocks for theories of highly correlated sys- 
tems. The recent success in incorporating such theories 
H into the basic apparatus used in electronic-structure 
calculations, means that exciting new progress in this 
area is close at hand. Here, we provide a surprisingly 
simple and elegant solution to the time-scale problem. 

A half a century ago, the framework and starting point 
for discussion of the time scale r for a localized spin in an 
electron sea was set by Korringa || , whose contribution 
has been promulgated at the textbook level for decades 
. In simple terms 1 /r gives the fractional rate at which 
a component of spin S representing a magnetic impurity 
(or a quantum dot) is changing due to the electrons in 
the conduction band (or in the leads of a quantum dot). 
It is given by an expression of the type 

- = ak B T, (1) 

T 

where a is a dimensionless constant, k B is Boltzmann's 
constant, and T is the temperature. 

The simplest model that contains the relevant physics 
is the Kondo model Hamiltonian X)"=i( efc i — ' ^0> 
where ki is a generalized quantum number for the i th 
conduction (or lead) electron, e^ is its energy, and 
is its spin operator. For a quantum dot, k is assumed 
to include the information on which lead is referred to. 
The exchange coupling J is taken to be independent of 
k, aside from the usual high-energy cutoff at energies fur- 
ther than D above or below the Fermi level. The Kondo 



temperature Tk is the single parameter characterizing 
the physics, and is of order {D/k B ) exp(— 1/Jp), where 
the density of states p is assumed constant. 

The factor of k B T in Eq. ([!]) arises because a change 
in S involves the absorption of a spin one electron-hole 
pair of zero energy. The phase space for such an object 
is doubly restricted by the Pauli principle, and produces 
the factor fdef(e)(l — /(e)) = k B T where /(e) is the 
Fermi function. Although Korringa's original derivation 
applied to a nuclear spin where the interaction with elec- 
trons is dipolar, it has been widely applied to impurity 
electron spins as well, using the Kondo model for the in- 
teraction. To lowest order in J, the quantity a in Eq. (|l|) 
is given by [| 

a = 7r(Jp) 2 . (2) 

One should note that (Q) represents a difference between 
rates and cannot be calculated directly by a second-order 
perturbation calculation. It rather involves the solution 
of a master equation, or, in many-body terminology, a 
sum of vertex corrections. For S — \, as would be appro- 
priate for a quantum dot, the vertex corrections change 
the result only by a factor of order unity. 

The application of these ideas to a biased quantum dot 
allows the phase space factor to be tuned in a continuous 
fashion by varying the bias V across the two leads of the 
dot, which for simplicity we assume to be symmetric un- 
der lead interchange. The phase space restriction factor 
in this case is given by \ v J de fi(e)(l - fv(e)), where 
the indices / and V designate which lead is referred to. 
For example, the Fermi functions /i(e) and /2(e) have 
Fermi levels displaced by ±eV72, respectively, where e 
is the magnitude of the electronic charge and V is the 
bias. The integral above can be evaluated analytically, 
with the result that Eq. (|l|) should be replaced by 

where 

F(x) = -(l + l + — - — + — 1 . (4) 

V ; 4 V 1 - e-* 1 - e~ x J V ' 

In writing Eq. (Q) we have sacrificed conciseness to fa- 
cilitate clarification of the origin of the terms in paren- 
theses, in terms of the wave function of the annihilated 
particle-hole pair. The first two terms arise when the two 



components of this object are in the same lead; in this 
case the existence of V has no effect on the result [||; 
the phase space for these processes is still constricted, 
and the contribution to 1/r is still small. For the third 
term, the particle is on lead 1 and the hole on lead 2; 
here the phase space is opened wide by V. Finally for 
the fourth term, the particle is on lead 2 and the hole on 
lead 1; the phase space is, aside from an exponentially 
small tail, closed off entirely as V is increased. So the es- 
sential physical feature deriving from Eq. (|4|) is that the 
factor of k B T in Eq. (|l|) is replaced at large V by \eV. 
The notion of an expanded phase space is implicit in the 
Anderson model calculation of r by Wingreen and Meir 
fiof and in a different context in the work of Kaminski et 
al. [ fill Here we have expanded the notion and expressed 
it more succinctly, because it is one of the cornerstones 
of what follows. 

We now turn to the calculation of a when there is a 
significant Kondo effect. The quantum dot model will 
also provide a framework for a precise definition of r. A 
primary feature of the Kondo effect is the appearance of 
a resonance in the the imaginary part of the transition 
or t matrix, which is centered at the Fermi level, and 
which extends logarithmically over many energy scales, 
the smallest of which is characterized by the Kondo en- 
ergy fee 7k- The transition matrix then determines the 
conductance G = I/V, where here I is the steady-state 
current induced by a constant bias voltage V. Thus a 
current measurement through a biased quantum dot pro- 
vides a direct measurement of the changes in the Kondo 
resonance induced by a time dependent J in the Kondo 
Hamiltonian. Specifically, we imagine that J, initially 
zero, is suddenly switched to its final value. The current 
I(t) will increase, rapidly at first, but for large times t 
should approach its final value of GV exponentially with 
a constant rate. The exponent for this final decay should 
be independent of the method of excitation, and hence 
the true characteristic formation rate of the Kondo reso- 
nance. 

Although the physics implicit in Eq. (||) and the order 
of magnitude it gives, has been used by several authors 
|^0| , ^T| to obtain time scales of the quantum dot, it is un- 
doubtedly incorrect to expect that this expression, valid 
to order J 2 only, can capture the time scale for a phe- 
nomenon like the Kondo resonance, whose defining terms 
first appear in order J 3 . Fortunately one can use previ- 
ous analytic results [fL2| to extract this fundamental time. 
What was found there, was that the current response to 
a stepped turning on of J would be the same as the equi- 
librium response to a time-dependent effective tempera- 
ture T c given by T c = T cot\vTrTt/2h. We start with the 
tautology SI /I = [dQnG)/dQnT e )] ST C /T C . Defining the 
fraction / as the finite difference 



and 



dQnG) 
lQ - d(lnT) : 



(6) 



we have approximately / = fo(cothirTt/2h — l)/2. As 
t becomes large, this approximation becomes exact, and 
we obtain 

/-►/oexpf-^Y (7) 

where the value of a defining r via Eq. (Q) is given by 

a = 7T, (8) 

instead of by Eq. (Q). When the phase space expansion 
factor (^) for finite V is restored 13 , then the predicted 
t should be given by Eq. (||). 

Our result for r, although analytically derived only in 
perturbation theory, contains no non-universal parame- 
ters like J or D, and so might be expected be valid in 
parameter ranges exceeding those implied in its deriva- 
tion. However, the expressions do not contain 7k, so 
corrections may appear when k B T and eV become much 
smaller than fee 2k- For example, one could suspect that 
a particle-hole pair would require an energy of ^fceTk 
in order to excite the singlet ground Kondo state, so 
that 1/r would vanish exponentially as ~exp(— Tr/T) 
for eV < k B T < k B T K instead of linearly with T. 
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(5) 



FIG. 1. Instantaneous current I through the quantum 
dot versus the time t after J was turned on. The units 
J u = (e 2 /nh)V and £k = ft/fcaTk- The solid, dashed, dotted, 
and dash-dotted curves correspond to the respective eV/ksTK 
values of 1, 5, 20, and 80. The temperature T = 1.5Tk. 
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We have made numerical simulations |14[] to test the 
response of the current to sudden increases in J at var- 
ious biases and temperatures. The essential feature is 
that J is suddenly increased from a negligible value to 
a substantial value, in this case given by Jp = 1/2tt. 
We found results for two temperatures: T — 1.5Tk and 
T = 15Tk, where Tk corresponds to the final value of J 
above. A sampling of the data is shown in Fig. 1. 

One sees clearly that the current versus time curves 
have the general features discussed above, with the time 
required for saturation decreasing with increasing bias V. 
For eV somewhat greater than fceTk, small, damped os- 
cillations occur jl5^[l7) at the non-trivial frequency given 
by Tiuj = eV. They are caused by transitions between 
the split Kondo peaks (SKP), which occur [flbjpl for 
eV greater than several times ^b^k- Such SKP oscil- 
lations should not occur at all for the completely sym- 
metric mechanism of suddenly turning on J. However, as 
described below, in order to make our simulations numer- 
ically feasible, it was necessary to include a small (asym- 
metric) mixed valent admixture, which is responsible for 
the excitation of the SKP oscillations in this case. 
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FIG. 2. Comparison of Eq. (g) with the simulation re- 
sults. The ordinate shows the rate l/tj/io as denned in the 
text in units of l/ix = /cbTk/Ti. For T = 1.57k, the dashed 
and solid lines show the theoretical curves for ax/io values 
of 1.6 and 2.0 respectively, while the squares show the simu- 
lation results. The dotted and dash-dotted lines and circles 
respectively show the corresponding results for T = 157k . 

The scheme we adopted for analyzing the data of 
Fig. 1, comprised first of determining the time, which we 
call tf, corresponding to the current fraction / defined 
earlier. We choose this fraction to be 1/10, an ampli- 



tude that is sufficiently large that the SKP oscillations 
do not spoil the analysis, and sufficiently small that the 
ultimate time scale has been reached approximately. One 
finds, by solving Eq. ([?]) for t, that tf is given simply by 
replacing r by tf in Eq. (|^) , provided that a is replaced 
by ctf = a/ln(/o//)- The amplitude /q was calculated 
directly both from the large-time asymptotes of the sim- 
ulations and from Eq. (|^), with similar results; it was 
found to be remarkably slowly varying with temperature 
in the temperature range of the simulations. Jl9| Assum- 
ing that a = 7r, we found taking / = 1/10 that a.f = 1.6 
for both T = 1.5T K and T = 15T K . In Fig. 2 we show 
our simulation results along with the theoretical predic- 
tion from Eq. (^) with ai/ w — 1.6. Better fits to the 
simulation can be obtained by taking ax/io — 1-8 (not 
shown) or aq/io = 2.0; these correspond to a values of 
1.1-7T and 1.27T respectively. Such a 10-20% difference 
from the prediction of Eq. (||) is consistent with the ex- 
pected accuracy of the simulation method, exacerbated 
by the fact that the simulations contained a small mixed 
valent contribution not present in the Kondo Hamilto- 
nian. The perturbative value of a — 0.03-7T predicted by 
Eq. (||) is definitely ruled out. The substantial agree- 
ment gives added support to our time scale predictions 
of Eqs. (0) through (||). The challenge of extending our 
results to the range k^T, eV <C fceTk still remains. 
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Note added on 27 September 2001. Since this work 
was originally posted, we became aware of work by 
A. Rosch, J. Kroha, and P. Wolfle, Phys. Rev. Lett. 
87, 156 802 (2001). These authors calculate the imag- 
inary part of the pseudo-fermion self-energy in time- 
independent NCA. If this quantity is closely related to 
the time-dependent definition, Eq. (7), of the Kondo 
state's formation rate, then their work would include cor- 
rections that go beyond the simple argument that led to 
Eq. (3), but which would automatically be included in 
our NCA simulations. It is possible that the small de- 
viation of the squares in Fig. 2 above the solid line for 
the highest three points could be due to this. However, 
the accuracy achievable by our data analysis is limited 
in this region and drawing any physical conclusion from 
this deviation might be straining these limits. Because 
the relative effect of the mixed valent contribution in- 
creases with increasing V, we are unable approach the 
asymptotic limit exemplified by Eq. (4) of the work of 
Rosch et al. We thank Achim Rosch for correspondence. 



3 



[1] D. Goldhaber-Gordon, et ai, Nature 391, 156 (1998). 
[2] S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 

hoven, Science 281, 540 (1998). 
[3] V. Madhaven et ai, Science 280, 567 (1998). 
[4] J. Li, W. D. Schneider, R. Berndt, and B. Delley, Phys. 

Rev. Lett. 80, 2893 (1998). 
[5] S. V. Savrasov, G. Kotliar, and E. Abrahams, Nature 

410, 793 (2001). 
[6] J. Korringa, Physica 16, 601 (1950). 
[7] C. P. Slichter, Principles of Magnetic Resonance 

(Springer, Berlin, 1990). 
[8] For example, see D. C. Langreth and J. W. Wilkins, 

Phys. Rev. B 6, 3189 (1972), Eq. (1.16). 
[9] A possible consequence of this fact has recently been ex- 
plored by P. Coleman, C. Hooley, and O. Parcollet, Phys. 

Rev. Lett. 86, 4088 (2001). 
[10] N. S. Wingreen and Y. Meir, Phys. Rev. B 49, 11040 

(1994). 

[11] A. Kaminski, Y. V. Nazarov, and L. I. Glazman, Phys. 
Rev. Lett. 83, 384 (1999). 

[12] P. Nordlander, et ai, Phys. Rev. Lett. 83, 808 (1999). 

[13] The implicit assumption that a itself is insensitive to V 
is consistent with its T independence, but should be a 
subject of future investigation. 

[14] The simulations used the noncrossing approximation, a 
reliable calculational scheme, whose time-dependent ver- 
sion we have adequately described and referenced pre- 
viously in Refs. Jl2 , jl5| . It is valid in the temperature 
range of our simulations, but breaks down for T <C Tk- 
In the same notation as we used previously in Ref. JlBj l, 
we took Anderson model parameters edot = — 2Pdot, with 
a parabolic band of width 9rdot- These give an equiva- 
lent Kondo model Jp = 1/27T and fceTk = 10 _3 Pdot. 
The equivalence between the models is good, but not 
perfect, because for computational-resource reasons, the 
small mixed valent overlap could not be reduced to zero 
by choosing a much larger |edot|- F° r times t < 0, the cou- 
pling J was for practical purposes turned off by taking 

Edot = — 5Pdot- 

[15] M. Plihal, D. C. Langreth, and P. Nordlander, Phys. 
Rev. B 61, R13 341 (2000). 



[16] P. Coleman, et ai, Phys. Rev. B (submitted); cond- 



matt/0108001 



[17] A. Schiller and S. Hershfield, Phys. Rev. B 62, R16 271 
(2000). 

[18] N. Sivan and N. S. Wingreen, Phys. Rev. B 54, 11622 
(1996). 

[19] Interestingly, fo can be evaluated analytically for the 
high-temperature region InT /7k S> 1, although such 
temperatures would be way out of the range used 
in our simulations. If we note that very early work 
[A. A. Abrikosov, Physics 2, 5 (1965)], valid in the 
high T range defined above, implies that (h/2e 2 )G — > 
37r 2 /161n 2 (T/T K ), then we find from Eq. || that f -> 
4/ln(T/Tk). More recent work [A. Kaminski, Y. V. 
Nazarov, and L. I. Glazman, Phys. Rev. B 62, 8154 
(2000)] suggests that V dependence might also be im- 
portant for comparably large values of V. 



4 



